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We carry out numerical simulation of the formation of U(l) global vortices in a first order phase 
transition in 2+1 dimensions in the presence of small explicit symmetry breaking. Bubbles of broken 
symmetry phase are randomly nucleated, which grow and coalesce. Vortices form at junctions of 
bubbles via standard Kibble mechanism as well as due to a new mechanism, recently proposed by 
us, where defect-antidefect pairs are produced due to field oscillations. In a simulation involving nu- 
cleation of 63 bubbles, with bias in phase distribution inside bubbles arising from explicit symmetry 
breaking, we find that not a single vortex/antivortex is produced via the Kibble mechanism, while 
the new mechanism leads to production of 104 vortices and antivortices. Even without biasing the 
phase distribution inside bubbles, the vortex production is completely dominated by this new mech- 
anism, which accounts for the production of about 80% of the vortices and antivortices, remaining 
20% being produced via the Kibble mechanism. We study the dependence of the effectiveness of 
the new mechanism on the magnitude of explicit symmetry breaking, as well as on the nucleation 
rate of bubbles. We also study the effect of damping on this mechanism and show that damping 
suppresses this mode of vortex production. 
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I. INTRODUCTION 

There has been a renewed interest recently in studying 
the formation of topological defects in phase transitions. 
It is clearly important to understand various processes 
which could be responsible for defect formation as then 
only one can hope to get a true understanding of the dis- 
tribution of defects and its time evolution. These issues 
are of great importance not only for condensed matter 
physics but also in the context of particle physics models 
of the early Universe where topological defects are sup- 
posed to play an important role in the evolution of the 
Universe and in the structure formation . 

In conventional studies only two distinct types of pro- 
cesses have been considered to be important in studying 
the formation of topological defects in phase transitions. 
One of these is based on thermal production of defects 
which leads to defect density which is suppressed by the 
Boltzmann factor ||^. The other process is based on the 
formation of a kind of domain structure during phase 
transition with defects forming at the junctions of these 
domains, and is typically known as the Kibble mecha- 
nism 1^. Kibble mechanism was originally proposed for 
studying defect formation in the early Universe. How- 
ever, the mechanism as such has complete general appli- 
cability and in fact has been recently verified by study- 



ingdefect formation in certain condensed matter systems 
[Though, there are nontrivial issues for the case 
when gauge fields are also present. See ref. Q and refer- 
ences therein.] 

An important aspect of the Kibble mechanism is that 
it does not crucially depend on the dynamical details of 
the phase transition. For example, the number of defects 
(per domain) produced via the Kibble mechanism de- 
pends only on the topology of the order parameter space 
and on spatial dimensions. Dynamics plays a role here 
only in determining the relevant correlation length, which 
in turn determines the domain size, affecting net number 
of defects produced in a given region. Still, the number 
of defects per domain is entirely independent of the dy- 
namics. [Apart from some special situations, e.g. in a 
very slow first order transition, see [Q .] 

Recently, a new mechanism for defect production has 
been proposed by us [P-p^ where the dynamics of the 
order parameter field plays a very important role. Here, 
defect-antidefect pairs are produced due to strong oscil- 
lations of the field. We showed that whenever the field 
passes through zero magnitude, while oscillating, (in a re- 
gion where the field is non- uniform), a defect-antidefect 
pair gets created. This essentially explains why the dy- 
namics plays a major role in this mechanism since the na- 
ture and the strength of field oscillations during a phase 



* E-mail : 
digal@iopb.stpbh.soft.net 
supratim@iopb.stpbh.soft.net 
ajit@iopb.stpbh.soft.net 



1 



transition will in general depend on various dynamical 
details, such as the order of the transition, time dura- 
tion of the quench etc. We find that the number of de- 
fects (say per domain) produced via this mechanism can 
vary strongly as different parameters of the transition are 
changed. 

This mechanism was first proposed by two of us jsj for 
systems in the presence of small explicit symmetry break- 
ing. Using numerical simulation of two bubble collisions 
in 2-1-1 dimensions for the case of U(l) global vortices, it 
was shown in ref. |^ that a very large number of vortex- 
antivortex pairs can be created via this mechanism due 
to field oscillations enhanced by the explicit symmetry 
breaking. For example, in one case it was found that 
ten, well separated, vortices and antivortices were cre- 
ated in a single two bubble collision. It was later shown 
by us in ref. that this mechanism was not restricted 
only to systems with explicit symmetry breaking, and is 
in fact very general. It applies even when there is no 
explicit symmetry breaking and is applicable to all sorts 
of topological defects. We also argued that this mecha- 
nism should be operative even in a second order phase 
transition involving quenching from a sufficiently high 
temperature. However, as this mechanism depends sen- 
sitively on dynamics, we found that defect production for 
the case of zero explicit symmetry breaking was not very 
prominent compared to the case with explicit symmetry 
breaking studied in ref. 1^. [At least for the range of 
bubble separations considered in ref. |8[|^.] 

Here we mention the work of Copeland and Safhn PJTTI 
who have considered the production of strings in gauge 
theories. They showed in ref. for Abelian Higgs 
model, that the geodesic rule may get violated in the 
collision of two bubbles and that vortex-antivortex pairs 
can form in the region of coalescence. In their case also 
the oscillations of field played crucial role, though due to 
the presence of gauge field the dynamics of the phase 6 
had extra features. For example, the presence of gauge 
fields provides a driving force for 9. The dynamics of 
vortex production in ref. |^ is, in this sense, similar to 
the case of explicit symmetry breaking discussed in ref. 

The studies in ref. were carried out by choosing 
specific initial bubble configurations. In order to study 
this mechanism in a realistic situation of a phase transi- 
tion, we carried out a full study of the eflFectiveness of this 
mechanism, for the case of zero explicit symmetry break- 
ing, by simulating a first order transition via random nu- 
cleation of bubbles . We showed there that for very 
low bubble nucleation rates (leading to large inter-bubble 
separation and hence very energetic collisions), the new 
mechanism becomes the dominant mode of production 
of well separated vortices. However, for large nucleation 
rates the vortex-antivortex pairs produced via this mech- 
anism consist of strongly overlapping field configurations 
which decay rapidly. Thus for large nucleation rates, 
number of surviving vortices will be roughly the same as 
expected from the Kibble mechanism. These results are 



important as they show that defect production in certain 
situations (i.e. very low nucleation rate) may be drasti- 
cally altered due to contributions of this mechanism. For 
example, this mechanism may completely dominate the 
production of cosmic strings, monopoles etc. in models of 
extended inflation in the early Universe, and may play an 
important role in the production of vortices in superfluid 
^He A - B transition 

From the results in ref. Q we expect that there is a 
special class of systems, those with small explicit sym- 
metry breaking, where this mechanism may be remark- 
ably effective in defect production during a phase transi- 
tion. There are many examples of systems with sponta- 
neously broken symmetries where the symmetry is also 
broken explicitly. In the context of particle physics, the 
Skyrmion picture of baryons in the context of chiral mod- 
els is an example where explicit symmetry breaking terms 
are needed to incorporate a non-zero pion mass. In fact 
it was this system where the role of explicit symmetry 
breaking in enhancing defect production was first dis- 
cussed by Kapusta and one of the author , though the 
underlying mechanism was different there. There is an- 
other example, in Particle physics, that of axionic strings 
arising from the breaking of the Pecci-Quinn symmetry, 
where explicit symmetry breaking is present [Q. It was 
suggested in ||^-[ic|] that this new mechanism may be es- 
pecially effective for the production of axionic strings due 
to the presence of explicit symmetry breaking. However, 
that argument does not seem to be correct since the scale 
of Pecci-Quinn symmetry breaking is very large (of order 
10^*^ GeV or so) compared to the scale of explicit sym- 
metry breaking, which arises at the QCD scale. Thus, at 
the Pecci-Quinn scale, explicit symmetry breaking can 
not play any role in the production of axionic strings, 
though this new mechanism may certainly contribute to 
their production. 

In condensed matter, nematic liquid crystals provide a 
simple example of such systems where the presence of ex- 
ternal electric or magnetic fields induces explicit symme- 
try breaking terms Pq | . Hence this new mechanism may 
play a dominant role there. Though there are subtleties 
in considering this system due to the fact that opposite 
orientations of the order parameter are identified. We 
will discuss this point later in Section V. 

From the point of view of these systems, and in general 
for systems with explicit symmetry breaking, it is impor- 
tant to know the actual contribution of this mechanism in 
a realistic phase transition. The study in ref. j|] consid- 
ered some specific bubble configurations and showed how 
a very large number of defects may be produced in a sin- 
gle two bubble collision. What one needs to know is the 
average number of defects produced via this mechanism 
when bubbles are randomly nucleated. We have already 
carried out such a study for the case of zero explicit sym- 
metry breaking, which is the case of widest possible appli- 
cations. [For example, the production of cosmic strings, 
magnetic monopoles etc., in the early Universe, and the 
production of vortices in superfluid '^He and *He in the 
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laboratory.] As the detailed dynamics, and hence the de- 
fect production via this mechanism, is very different for 
the systems with explicit symmetry breaking, we need 
to carry out detailed simulations of the phase transition 
for this case as well. Certainly, for baryon (Skyrmion) 
production in the heavy-ion collisions, as well as for the 
production of liquid crystal strin gs in the presence of ex- 
ternal fields, the results of ref. [|lO[ (for the zero sym- 
metry breaking case) are not much useful. We should 
expect that a much more dramatic role will be played by 
the new mechanism for these systems, even in the generic 
situation of a phase transition. 

There is one more reason why this mechanism may be 
most dominant in the situation with explicit symmetry 
breaking. If bubble nucleation is not happening at very 
large temperatures (compared to the explicit symmetry 
breaking scale) then the distribution of phases inside dif- 
ferent bubbles will get biased. This bias will suppress 
vortex production due to Kibble mechanism. For ex- 
ample, in our simulation with 63 bubbles (discussed in 
Section IV), essentially all the bubbles had 9 which was 
either less than 7r/2 or greater than 37r/2 (for the case 
when 9 = corresponds to the true vacuum). Due to 
this, not a single vortex was produced via the Kibble 
mechanism. Bias in 9 towards value zero does have ad- 
verse effect even on the new mechanism due to reduced 
potential energy which affects the flipping of the field. 
However, wall oscillations will be expected to be much 
larger now due to smaller 9 gradients between colliding 
bubbles. This should help vortex production via the new 
mechanism. What we find is that all the vortices and 
antivortices (a total of 104 in this case) were produced 
via the new mechanism. This is about the same average 
vortex production per bubble via the new mechanism as 
we obtain when 9 distribution inside bubbles is chosen to 
be uniform (though in that case there are some vortices 
produced via the Kibble mechanism as well). 

We also would like to know how this average defect pro- 
duction is affected by dynamical details, such as the mag- 
nitude of explicit symmetry breaking, nucleation rate, 
presence of damping etc. We carry out a complete study 
of all these issues in this paper. We find that defect pro- 
duction increases with increasing magnitude of explicit 
symmetry breaking. Defect production also increases for 
lower nucleation rates due to bubble collisions being more 
energetic. This is the same behavior which was observed 
in ref. for the case of zero explicit symmetry break- 
ing. Presence of damping suppresses field oscillations and 
hence results in lower number of defects. However, com- 
pared to the case of zero explicit symmetry breaking dis- 
cussed in ref. ||l^ where strong damping makes this mech- 
anism completely ineffective, for the present case we find 
that some pairs are always produced via this mechanism 
(as long as explicit symmetry breaking is not too small) . 
This results is very significant for the cases of experi- 
mental interest. For example, in the case of Skyrmion 
production in heavy- ion collisions, dissipative effects will 
always be present in the expanding plasma. This is es- 



pecially true for the case of liquid crystal strings where 
the dynamics is completely dominated by dissipation. 

The paper is organized in the following manner. The 
second section reviews the essential physical picture of 
the new mechanism and discusses the numerical tech- 
nique. Section HI discusses the results of the simulation 
of the phase transition by random nucleation of bubbles 
with unbiased 9 distribution inside bubbles. In Section 
IV we present results for simulation with 9 distribution 
inside bubbles being biased according to the explicit sym- 
metry breaking term. The dependence of vortex pro- 
duction on nucleation rate and the magnitude of explicit 
symmetry breaking is discussed in Section V. The effect 
of damping on this mechanism is discussed in Section VI 
and conclusions are presented in Section VII. 

II. PHYSICAL PICTURE OF THE MECHANISM 
AND NUMERICAL TECHNIQUE 

We first briefly recall the essential features of this 
mechanism Q. Consider the case of spontaneous break- 
ing of a global U(l) symmetry, with the order param- 
eter being the vacuum expectation value of a complex 
scalar field <&. The vacuum manifold is a circle S^. After 
the phase transition, <i> will have a non-zero magnitude 
but its phase 9 may vary spatially. Consider a region 
of space where 9 has some small variation between two 
points A and B. Now suppose that the magnitude of the 
field undergoes strong oscillations in a small region be- 
tween points A and B. If $ passes through zero magni- 
tude while oscillating in a region, then it is easy to see 
that 9 in that region will discontinuously change to 9 + tt. 
We call it the flipping of <&. By drawing the distribution 
of 9, with flipped orientation in the middle portion be- 
tween A and B, one can see that a vortex and antivortex 
pair has formed in the opposite sides of the line joining 
A and B, see ref. for details. These considerations can 
easily be generalized to other defects |9j. 

Numerical techniques used in this paper are the same 
as used in the earlier papers |§-|lO|. We study 2+1 dimen- 
sional case, with the Lagrangian taken as the following. 

This Lagrangian is written in terms of appropriately 
scaled coordinates, and a dimensionless complex scalar 
field $, with (f> and 9 being the magnitude and the phase 
of Value of e is taken to be 0.05. This Lagrangian 
describes a theory where U(l) global symmetry is spon- 
taneously broken, except for the presence of last term 
which breaks U(l) explicitly. 

For the case of k = 0, the process of vortex creation via 
bubble nucleation has been studied in detail in ref. jlj] 
for the Kibble mechanism, and in ref. |||,|l^ for the new 
mechanism. At zero temperature, the phase transition 
takes place via nucleation of bubbles of true vacuum in 
the background of false vacuum via quantum tunneling. 
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The bubble profile 4> is obtained by solving the Euclidean 
field equation 

^ + ^^-m=0 (2) 
dr^ r dr 

where 1^(0) is the effective potential in Eq.(l) (with 
K = 0) and r is the radial coordinate in the Euclidean 
space. In the Minkowski space initial profile for the bub- 
ble is obtained by putting t — Q va. the solution of the 
above equation. 9 takes a constant value inside a given 
bubble. 

The bubble profile found by solving Eucledian equation 
of motion Eq.(2) for k = provides an adequate starting 
bubble configuration even for small non-zero values of k 
as we are only interested in vortex formation which hap- 
pens after bubbles expand and coalesce. This was the 
approach used in our earlier work |^ for the case with 
K 7^ 0. However, there is a nontrivial issue here. [This 
was pointed out to us by Copeland and Saffin] . The stan- 
dard method of finding the bounce solution for the 
present case with /t 7^ 0, works fine if the phase 9 takes 
value or TT inside the bubble. For any other choice 
of 0, the method does not work with the conventional 
approach. It appears that if at all there is any bounce 
solution for arbitrary 0, it must require a non-uniform 9 
inside the bubble 

For the purpose of the present paper, we will ignore 
this issue. One approach can be that if one was consid- 
ering thermal nucleation of bubbles, then if the explicit 
symmetry breaking term is very small compared to the 
scale set by the temperature, then one can ignore it while 
considering the nucleation process itself. Thus we will 
continue to use the bubble profile obtained for the case 
of K = and evolve the configuration with non-zero value 
of K. The profile of the bubble calculated with n — has 
the asymptotic value </> = which is the local minimum of 
the effective potential. Hence it is suitable for nucleating 
the bubbles in the background of false vacuum, = 0. 
This suggests that the explicit symmetry breaking term 
should be chosen such that the false vacuum remains the 
same (f) — so that there is no mismatch between the 
false vacuum and the asymptotic value of $ of differ- 
ent bubbles. The specific form of the explicit symmetry 
breaking term we use is motivated by this consideration 
as well as by simplicity. It is important to mention that 
our results are not sensitive to the form of the explicit 
symmetry breaking term, as long as the true vacuum 
(with non-zero value of k) is non-degenerate. We have 
verified that similar enhancement in vortex production 
results with other types of symmetry breaking terms as 
well, for example with K(f)cos9. 

Even though we are neglecting the effects of explicit 
symmetry breaking in solving Eqn.(2) for the bubble pro- 
file, its effects should be taken into account when 9 in- 
side bubbles is being determined [9 inside a given bubble 
taken to be uniform). For k — one takes all values of 
9 to be equally probable. However, with k 7^ one will 



expect a bias in 9 distribution inside bubbles as bubbles 
with smaller values of 9 will have smaller free energy and 
hence larger nucleation probability. We take into account 
this bias in 9 distribution inside bubbles later in Section 
IV. We show there that, in a simulation involving nucle- 
ation of 63 bubbles, bias in 9 distribution leads to all bub- 
bles having 9 configurations such that not a single Kibble 
vortex/antivortex is produced. At the same time vortex 
production due to the new mechanism is essentially un- 
affected by the biased 9 distribution. However, in order 
to have a conservative estimate of the importance of the 
new mechanism compared to the Kibble mechanism, we 
first study situation where bubbles are nucleated with 
uniform probability of 9 for which Kibble mechanism con- 
tribution is well understood. This approximation will be 
consistent when the tilt is small compared to the temper- 
ature scale at which bubble nucleation is supposed to be 
happening. [Though, in that situation one should worry 
about validity of classical evolution using field equations 
as for large temperatures vortex-antivortex pairs can be 
thermally nucleated. We will consider this point later 
in Section IV. For now we will continue using uniform 
probability distribution for 9 inside bubbles.] 

Rest of the numerical technique is the same as used 
in earlier works [p|-p^. The field evolution is carried out 
using field equations obtained from Eqn.(l). The simu- 
lation of the phase transition is carried out by nucleat- 
ing bubbles on a square lattice with periodic boundary 
conditions, i.e on a torus. Lattice spacing in spatial di- 
rections was taken to be Aa; = 0.16, with time step At 
equal to /\x/^/2. Simulations were carried out on an HP- 
735 workstation, a Silicon Graphics Indigo 2 workstation, 
and an HPK-260 workstation at the Institute of Physics, 
Bhubaneswar. 

When evolved with these equations, bubbles grow, co- 
alesce, and vortices are formed in the intersection region 
of three or more bubbles by Kibble Mechanism, as well 
as by the decay of bubble walls via the new mechanism. 
Note that, when k is non zero then 6* = is energetically 
preferred (for our choice of explicit symmetry breaking 
term). So when a bubble with non-zero 9 is nucleated, 9 
inside the bubble evolves towards 9 — and eventually 
starts oscillating about 6* = with decreasing amplitude. 
This non-trivial dynamics of 9 is one feature which makes 
the present case qualitatively different from the case of 
K = considered in ref. [Q. We explain this in some 
detail in the following by considering a two bubble colli- 
sion. 

When two bubbles collide there is always (f) oscillations 
in the coalesced portion of the two bubbles. Magnitude 
of this oscillation depends on the 9 difference between the 
two bubbles as well as on the spatial separation between 
the two bubbles. As described in ref. [|l^, (j) oscillations 
are more prominent when 9 difference between the bub- 
bles is small. Large 9 difference leads to large gradient 
energy in the coalesced region which suppresses (p oscil- 
lations in that region. When ip oscillations have suffi- 
cient amplitude then $ can overshoot the value $ = 
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by climbing over the potential hill. Clearly, this leads to 
a change in 9 in that region to + tt, i.e. to flipping of 
$. If this oscillation, and the associated flip, happens 
inside a region where is spatially varying, then a pair 
of vortex-antivortex forms in this single flip (as explained 
at the beginning of this section). Note that so far in this 
picture we did not use anywhere the presence of explicit 
symmetry breaking (or, the rolling of 9 to zero). One will 
expect that every time this type of situation arises (even 
in the course of a single two bubble collision) , it will lead 
to a vortex-antivortex pair getting nucleated. In fact, 
this is where the explicit symmetry breaking term in the 
effective potential plays a crucial role, as we describe be- 
low. 

Consider two bubbles with associated values of 9 be- 
ing TT + a and tt — a respectively, with a < tt/2. When 
these two bubbles collide, according to geodesic rule, 9 
in the coalesced portion takes the value 9 = tt (for this 
a should be small as 9 rotates towards zero). Since 9 is 
rolling towards zero in both the bubbles we always en- 
counter = TT as we go from one bubble to the other 
and there is a gradient which keeps increasing due to 
rolling of 9. It helps to think in terms of an arc spreading 
around the value tt in the plot of the effective potential 
(see Fig.l) which represents the variation of 9 from the 
center of one bubble to the center of second bubble. Both 
ends of this arc (p and q) will roll down towards zero, and 
as the arc is constrained to go through the value tt (due 
to continuity), it will lead to larger 9 gradients. Now if 
$ flips in this region due to oscillation then we will get 
a vortex-antivortex pair. The vortex and the antivortex 
in this pair move away from each other as 9 is zero in 
between them while 9 in the outer directions is tt hav- 
ing larger potential energy. This is how the first pair is 
created. 



oscillation around zero in the bubbles. [9 in each bubble, 
while rolling down towards 0, overshoots leading to non 
zero 9 gradient in that region.] If 4> oscillations are still 
strong enough so that they lead to yet another flip of 
$ then another vortex-antivortex pair will be nucleated. 
This process of successive pair creations continues till we 
have 9 oscillations maintaining 9 gradient in the coalesced 
portion and large (j) oscillations flipping Creation of 
multiple pairs in this manner was demonstrated in ref. 

, where it was shown that as many as ten vortices and 
antivortices were created in a single, two bubble collision. 

We have emphasized that if there is no flipping of <& 
then there will be no pair creation of vortices. Also, the 
amplitude of (j) oscillation in the coalesced region depends 
on the separation of bubbles, with larger separation lead- 
ing to larger amplitude. This may lead one to think that 
if the bubbles are at the closest separation then there will 
be no large oscillation of (j) so there may not be any pair 
creation at all. This was indeed true for the case of k = 
discussed in ref. However, this is not true for the 

present case, with k ^ 0. Consider for example a configu- 
ration which leads to = tt in the coalesced region. The 
gradient energy density in 9 variation in the coalesced 
region will keep increasing as 9 in both the bubbles rolls 
towards 0. This, implies large energy in the coalesced 
portion due to explicit symmetry breaking term which 
can drive (j) to climb the potential hill and lead to the fiip- 
ping of <&. [At least when k is not too small, see below.] 
Therefore, there is at least one pair of vortex-antivortex 
pair nucleated for appropriate initial 9 even for smallest 
separation of the bubbles. As we will see later, even in 
the presence of strong damping one pair is still nucleated 
(in contrast to the case of ref. [Q), though successive nu- 
cleation may be suppressed depending on the magnitude 
of damping. 




Fig.l 



FIG. 1. Plot of the effective potential. 9 variation in the 
physical space is taken to correspond to the solid arc. The 
end points, p and q, of the arc roll towards 9 = 0, but the arc 
is constrained to go through 9 = tt due to continuity. 

We note now that even though 6* is in the coalesced 
portion, there is still a 9 gradient which arises due to 9 



III. SIMULATION OF THE PHASE TRANSITION 

The preceding section explains how vortex-antivortex 
pairs are created in a two bubble collision in the presence 
of explicit symmetry breaking. It was shown in ref. |^ 
that a large number of vortex-antivortex pairs can form 
in a two bubble collision for appropriately chosen initial 
field configuration of bubbles. Clearly, in an actual phase 
transition, such large enhancement may not be expected 
due to randomness in values of 9 inside the bubbles as 
well as in their separations. This is especially so because 
the vortex production via this mechanism depends on the 
dynamical details of the bubble collisions and depends on 
parameters such as the magnitude of symmetry breaking 
term in the effective potential. Therefore, to find actual 
enhancement in the vortex production in this case one 
has to carry out the full simulation of the phase transi- 
tion. We describe the results of such simulations in this 
section. 

In this section we will present results with the value of 
K — 0.015. Here, we will neglect the explicit symmetry 
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breaking term in determining 9 inside nucleated bubbles. 
In the next section we will take into account the effect of 
explicit symmetry breaking in biasing the 6 distribution. 
We will see that the average defect production per bub- 
ble, via the new mechanism, remains almost unaffected, 
while defect production via the Kibble mechanism be- 
comes strongly suppressed for biased 9 distribution case. 

Nucleation rate was suitably chosen so that there were 
all together seven bubbles nucleated in the whole lattice. 
Physical size of the lattice was taken to be 192 x 192. 
As we mentioned earlier, the location of bubbles in the 
lattice as well as the values of 9 inside them were chosen 
randomly. Number of bubbles was chosen to be small (by 
taking small nucleation rate) so that bubble collisions are 
energetic which helps in vortex-antivortex pair creation. 
Also the bubbles are nucleated only in a short time dura- 
tion so they are roughly of same size when they collide. 
These bubbles then expand and coalesce. 

Fig. 2a shows the plot of at an early stage when bub- 
bles have just started coalescing. [In all the plots of $, 
the orientation of an arrow from the positive x axis cor- 
responds to 9 and the length of the arrow is proportional 
to (j).] Rotation of 9 can be clearly seen inside the bub- 
bles. Due to (j)^ dependence in the k term, 9 rotation is 
slowest near the bubble walls. Note that 9 in the central 
region of many bubbles has almost rotated to zero, the 
only memory of the initial values of 9 remaining near the 
walls. Interestingly, distribution of 9 at walls of different 
bubbles plays the most dominant role in the production 
of all the vortices observed in the simulations. 



The location of the vortices was determined by using 
an algorithm to locate the winding number. As the phase 
transition nears completion via the coalescence of bub- 
bles, magnitude of $ becomes non-zero in most of the re- 
gion with well defined phase 9. We divide each plaquette 
in terms of two (right angle) triangles and check, for each 
such triangle, whether a non-zero winding is enclosed. 
For this purpose we use the geodesic rule to determine 
configuration in between two adjacent lattice points, 
see ref. jl^ for details. Windings are detected only in 
regions where the magnitude of is not too small in a 
small neighborhood of the triangle under consideration. 
If $ is too close to zero in a region then that region is still 
mostly in the false vacuum and there is no stability for 
any windings present there. After getting probable loca- 
tions of vortices using the above algorithm, we check the 
region containing each vortex/antivortex, using detailed 
phase plots and surface plots of (j) to check the winding of 
the vortex, and select only those vortices which have well 
defined structure. By checking similar plots at earlier as 
well as later time steps we determine whether the vortex 
was produced due to oscillation, and subsequent flipping, 
of $, or via the Kibble mechanism. 
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FIG. 2. (a) Plot of $ at an early stage when bubbles have 
just started coalescing. In all the plots of $, the orientation 
of an arrow from the positive x axis corresponds to and 
the length of the arrow is proportional to (f>. (b) $ plot in 
a region where there is no winding present at this stage, (c) 
An elongated region towards left of the center is oscillating, 
and a pair has been created, (d) Another pair has been cre- 
ated in this region due to continued oscillations, (e) $ plot 
of the entire lattice at an stage when maximum number of 
vortices and antivortices are present. Altogether, there are 
four vortex-antivortex pairs present at this stage. 

Figs. 2b, 2c and 2d show a region where two vortex- 
antivortex pairs form due to field osciUations. Fig. 2b 
shows the plot of $ at t = 56.57. There is no winding 
present in the region at this stage. Fig.2c shows the plot 
at t = 62.22. Wc sec that an elongated region towards 
left of the center is oscillating, and a pair has been cre- 
ated. These oscillations continue, leading to creation of 
another pair, as shown in Fig. 2d. Fig. 2c shows the <I> 
plot of the entire lattice at this stage which corresponds 
to the situation when maximum number of vortices and 
antivortices are present. There are four vortex-antivortex 
pairs present at this stage. One of the pairs is near x = 
36, y = 96, second pair is near x = 36, y = 132, third 
pair is near x = 136, y = 72, and fourth pair is near x = 
116, y = 144. 

In this particular simulation a total of (time inte- 
grated) 14 vortices and antivortices were formed. Out 
of these, only two are formed due to Kibble mechanism. 
This is consistent with the expected number density of 
1 /4 per bubble (for U(l) vortices in 2 spatial dimensions). 
Remaining twelve vortices and antivortices are formed 



7 



due to this new mechanism of coalesced bubble wall de- 
caying by pair creation. This gives the average number 
of defects per bubble via the new mechanism to be about 
1.7. Thus in this case this mechanism is completely dom- 
inant over the Kibble mechanism. Average number of 
all defects per bubble in this case is then 2.0 which is 
eight times larger than the number expected if vortices 
were only produced due to the Kibble mechanism. In 
all our simulations we have found strong enhancement of 
the number density of defects per bubble in the presence 
of explicit symmetry breaking. We should mention that 
apart from these 14 vortices and antivortices there were 
four additional pairs (i.e. eight vortices and antivortices) 
formed where vortices and antivortices were not well sep- 
arated. Still windings of 9 corresponding to these pairs 
could be clearly seen in plots. We do not count such 
pairs as these defects annihilate quickly. Such strongly 
overlapping pairs were found in all simulations. 

There were some intriguing features observed in these 
simulations which we describe in the following. Fig. 3a 
shows the field configuration at t = 42.42, where we 
see a well separated vortex-antivortex pair. The vortex 
and the antivortex move towards each other, as shown 
in Fig. 3b. Due to asymmetric field configuration, the 
motion of the vortex-antivortex pair is not along the di- 
rection of their separation, rather the pair moves upward 
as a whole. The pair has annihilated hy t — 50.91, as 
shown by the absence of any windings in Fig. 3c. Subse- 
quent plot aX t = 53.74 in Fig. 3d shows that the vortex- 
antivortex pair has been re-created. Interestingly, how- 
ever, this time the location of the vortex and the antivor- 
tex are reversed. It is as if the vortex and the antivortex 
have passed through each other. Also, by checking sub- 
sequent plots, we find that the motion of the pair as a 
whole is now downward. This new pair annihilates later. 
This pair, thus, represents a bound state of a vortex- 
antivortex system. There were several vortex-antivortex 
pairs for which this cycle of annihilation and re-creation 
was seen. In all the cases, the vortex and the antivortex 
were exchanged after first annihilation and re-creation. 
For some cases, the cycle was repeated, but when the 
pair was re-created second time, the vortex and the an- 
tivortex scattered back, instead of going through each 
other. 

This process of the vortex and the antivortex pass- 
ing through each other in the first cycle of annihilation 
and re-creation should be contrasted with the results in 
ref. |lj,|2^ for vortex-antivortex scattering. For global 
defects it was argued in fl^ that for vortex-antivortex 
collisions which are dominated by gradient energy con- 
siderations, the vortex and the antivortex should retrace 
their original paths (i.e. they should bounce back) in a 
process of annihilation in a head-on collision and subse- 
quent re-creation. Numerical simulation of gauged vor- 
tices in ||2^ showed that in such collisions the vortex and 
the antivortex bounce back as long as their kinetic ener- 
gies are not very large, but they go through each other 
for extremely energetic collisions. In light of these stud- 



ies, our numerical results of passing through of vortex 
and antivortex suggests that the collision is sufficiently 
energetic that gradient energy considerations are not rel- 
evant any more, like the results in ref. |20| for large ki- 
netic energies. An estimate of the velocities from Figs. 3a 
and 3b shows the velocity of the vortex to be between 

0. 9 and 1 and that of the antivortex, between 0.8 and 

1. The uncertainty in the velocity arises from the un- 
certainty in the locations of the centers of the vortices. 
These are indeed very large velocities. One would expect 
such large velocities to arise here due to the domain wall 
(arising due to the concentration of 9 gradient) stretching 
between the vortex and the antivortex. Scattering back 
of vortices in subsequent annihilations and re-creations 
is then understood as due to velocities of vortices being 
not too large, so the arguments of ref. should be ap- 
plicable. This scattering back of vortex-antivortex thus 
provides an illustration of the physical arguments given 
in ref. [l9|] . We mention here that vortex-antivortex an- 
nihilations were seen in k = case also [|l^, but the 
vortex-antivortex pair was never re-created there. Pre- 
sumably, the velocities at annihilation are much larger 
for the present case of k 7^ 0. 
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Fig.3d 

FIG. 3. (a) A well separated vortex-antivortex pair, (b) 
The vortex and the antivortex are moving towards each other, 
while the pair, as a whole, moves upwards, (c) Pair has anni- 
hilated by this time, as shown by the absence of any winding 
in this plot, (d) Vortex-antivortex pair has been re-created. 
Note that the positions of the vortex and the antivortex have 
been exchanged. 

Another interesting result we find is for the cases of cer- 
tain two bubble collisions which lead to vortex-antivortex 
pairs which are not well separated. However, due to pres- 
ence of other bubbles, density waves arising from other 
bubble collisions propagate through the region where a 
given pair is getting formed. This density wave then 
strongly affects the (/> oscillations in that region and effec- 
tively separates the vortex and the antivortex. Thus, due 
to the effects of density wave, a pair which was strongly 
overlapping and was going to annihilate soon, ends up 
in well separated and well formed vortex and antivortex. 
This again shows that this new mechanism of pair cre- 
ation strongly depends on the dynamical details of the 
transition. A density wave which as such can not con- 
tribute to defect production, leads to enhancement in the 
production of well separated defects. Similar effects were 
also found in ref. llToi. 
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IV. DEFECT PRODUCTION WITH 
BIASED B DISTRIBUTION 



In this section we carry out a simulation as in Section 
III, but now with inclusion of a suppression factor in the 
nucleation rate for bubbles with a given d inside, result- 
ing from the explicit symmetry breaking. As we will see, 
Kibble mechanism vortices get strongly suppressed. To 
get an idea of this suppression, we carry out a larger 
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simulation with nucleation of 63 bubbles, while average 
inter-bubble separation is the same as in Section 111. This 
enables us to make direct comparison of average defect 
production per bubble for the two cases (as defect pro- 
duction via the new mechanism depends on average sep- 
aration of bubble nucleation sites). 

We estimate the suppression factor for bubble nucle- 
ation, with a given 9 inside, in the following manner. 
For this we assume that bubble nucleation is happening 
at some finite temperature T. (Even though the bubble 
profiles calculated in Section 11 correspond to T = case, 
we continue to use the same profiles. The reason for this 
is that we want to change only one factor here from the 
simulation of Section III, that is the 9 suppression fac- 
tor. In any case, these are thick wall bubbles for which 
the profiles for T = and T 7^ do not differ much. Also, 
actual profile of bubbles which is relevant for us is at the 
time when bubbles start coalescing. We are working with 
low nucleation rates, so bubbles expand by large amount 
before coalescing. Thus initial difference in the profiles 
becomes irrelevant.) This will lead to an additional fac- 
tor Tg in the nucleation rate of the bubbles which arises 
due to the dependence of the effective potential on 9. We 
have 

Te = e-^C')/^ (3) 

where F{9) is the contribution to the free energy of the 
bubble arising from the explicit symmetry breaking term. 
For T we use the constraint that T should be less than 
typical energy required to create a vortex-antivortex pair, 
otherwise use of field equations for the evolution will be- 
come suspect. We thus estimate the energy of a pair of 
vortex-antivortex at various separations. For this we use 
a code for minimization of energy as was used in ref . . 
To determine the energy of a vortex-antivortex pair, at 
a given separation, one starts with a trial profile for the 
pair. Field configuration is then fluctuated, while fixing 
the centers of the vortex and the antivortex, and energy 
is minimized. The configuration with the lowest value 
of energy is finally accepted as the correct profile of the 
vortex-antivortex pair and corresponding energy is taken 
as the energy of the pair. For details of the numerical 
technique, we refer the reader to ref. |2|]. 

We find that at separation of about (between the 
centers of the vortices) the energy of the pair, Epair, is 
about 4.2 ran, while at separation of 2m^ and 3m^^, 
Epair increases to about 5.8 rriH and 6.0 friH respectively. 
For the consistency of using field equations for the evo- 
lution, we require T < Epair- To get F{9), we integrate 
the explicit symmetry breaking term in Eqn.(l) for the 
profile of the bubble (determined from solving Eqn.(2)). 
We find F{9) ~ —3.6 mHCos9. We choose a sample 
value of T (satisfying the constraint that T < Epair) 
as T ~ 1.8 niH- Resulting value of F{9)/T is then sim- 
ply equal to —2.0 cos9. With this, we use the following 
(suitably normalized) expression for Tg for determining 
bias in 9 distribution in bubble nucleations. 



Te = e2(™''^-i) (4) 

With this choice of normalization, Yg is 1 for ^? = 
and is lowest, equal to e^^, for 9 — n. 

We then nucleate bubbles, with inclusion of this fac- 
tor for 9 suppression. To have reasonable statistics (due 
to strong suppression of Kibble vortices in this case) we 
carry out simulation over a larger physical region with 
size equal to 576 x 576. Total number of bubbles nucle- 
ated was 63 which lead to about the same average bubble 
separation as in Section III. Just as in Section 111, bub- 
ble locations are determined randomly, but now 9 inside 
bubbles is chosen with the weight factor Tg as give in 
Eqn.(4). 

Resulting 9 distribution of the bubbles is shown in 
Fig. 4. Solid dots show the frequencies P{9) of differ- 
ent values of 9 obtained in the simulation, normalized so 
that the maximum value of P{9) is 1.0. Solid curve is 
the plot of Tg. Most important thing to note in Fig. 4 
is an almost complete absence of 9 points in the range 
37r/2 > 9 > tt/2. In fact, there were three bubbles with 9 
being slightly less than 37r/2, or slightly more than tt/2. 
However, we verified that very quickly, during the evo- 
lution and much before any other bubble collides with 
these, 9 in these bubbles rotates (due to explicit symme- 
try breaking) and also falls in the range 37r/2 < 9 < tt/2. 
It is then immediately obvious that there is no possibil- 
ity of any Kibble mechanism vortex forming by coales- 
cence of these bubbles as no winding can be generated 
by following the geodesic rule. We have also verified it 
explicitly, by considering detailed field configurations of 
the vortices. 
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FIG. 4. 9 distribution inside bubbles. Dots show the nor- 
malized frequency P{9) of 6 inside different bubbles nucleated 
in the simulation. Solid curve is the plot of Tg. 

We get a total of 52 vortex-antivortex pairs (i.e. 104 
vortices and antivortices), all due to the new mechanism. 
This gives the average defect production per bubble to be 
about 1.65 which is essentially the same as the average 
number of defects per bubble via the new mechanism, as 
found in the simulation in Section 111. Thus biased 9 dis- 
tribution does not seem to adversely affect defect produc- 
tion via this new mechanism. This results is important 
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as it shows that in certain situations, with exphcit sym- 
metry breaking, this mechanism may not just dominate, 
it may be the only mechanism active for defect produc- 
tion. Effect of biased 6 distribution on this mechanism 
can be understood as follows. Suppression of larger val- 
ues of leads to smaller potential energy being stored 
in the collision regions, and this should adversely affect 
wall oscillation and hence defect production. However, 
now average 6 difference between coalescing bubbles will 
be smaller. As mentioned earlier, this helps in increasing 
wall oscillations, as less energy is spent in overcoming 9 
gradients. This should lead to increase in defect produc- 
tion. Final defect production will be a combination of 
both these factors. 

In next sections, we will study effects of changing nu- 
cleation rate, magnitude of explicit symmetry breaking, 
and presence of damping, on defect production via this 
new mechanism. As 9 suppression does not seem to af- 
fect this mechanism much, we will not include its effect 
in these sections. This also helps in making conservative 
comparison of this mechanism with the Kibble mecha- 
nism, when these parameters are varied, as in absence of 
any bias in 9 distribution Kibble vortices are also unsup- 
pressed. 

V. DEPENDENCE OF DEFECT PRODUCTION 
ON NUCLEATION RATE AND k 

The results discussed in last two sections show that 
there is a strong enhancement in the number of defects 
produced in a phase transition due to this new mecha- 
nism. In this section we will study the dependence of 
this enhancement in defect production on the magnitude 
of symmetry breaking coefhcient k and on the nucleation 
rate. We start by discussing the effects of changing the 
nucleation rate, which in effect changes the average sep- 
aration between bubbles. 

We use K = 0.015 as in Section III, and keep it fixed 
while we carry out simulations with different values of 
the nucleation probability. As we are now using uniform 
probability for 9 distribution inside bubbles, we will com- 
pare the results of this section with those in Section III. 
Bubble nucleation probability used for the case discussed 
in Section III lead to seven bubbles nucleated in the lat- 
tice. We have studied two other cases by changing the 
nucleation probability. In one case the net number of 
bubbles nucleated was 12 and in the other case number 
of bubbles was 20. 

For the case of 12 bubbles we found a total of 18 vor- 
tices. Out of these, 3 pairs of vortex-antivortex formed a 
cluster. For this simulation the total defect density per 
bubble comes out to be 1.5 which is less then the case 
of seven bubbles discussed in the previous section. This 
is consistent with the physical picture of this mechanism 
we discussed earlier where we mentioned that energetic 
bubble collisions will lead to larger oscillations of (f> in 
the coalesced portions of bubble walls and consequently 



to larger number of vortex-antivortex pairs getting pro- 
duced. Twelve bubble case leads to smaller value of av- 
erage separation between bubbles compared to the seven 
bubble case and hence less energy for (f) oscillations. 

As for the case of seven bubbles of Section III, 
here also we find several additional pairs of vortex- 
antivortex which correspond to strongly overlapping 
vortex-antivortex configurations. These did not separate 
and lasted only for a short time. There were three such 
pairs in the present case of twelve bubbles but we did not 
count them in the total number of vortices. 

For the second case we take the value of nucleation 
probability so that there are twenty bubbles nucleated in 
the whole lattice. Since the average separation between 
the bubbles is smallest in this case compared to the pre- 
vious two cases, the field oscillations flipping $ will be 
weakest in this case. Thus we expect that the pair pro- 
duction mechanism will not lead to many defects in this 
case. This is precisely what we find as the total number 
of vortices and antivortices in this case comes out to be 
14. This gives the number density of all defects per bub- 
ble to be equal to 0.7. Again, apart from these 14 defects, 
we found six vortex-antivortex pairs which are strongly 
overlapping and therefore are not counted in determin- 
ing the net number of defects. These results are similar 
to those found in ref. |jl^ where larger nucleation rate 
lead to a smaller number of well separated pairs being 
produced. 

To study the effect of k on the number density of de- 
fects we took the initial field configuration to be that 
corresponding to the case of twelve bubbles mentioned 
above. We then carry out simulations for various values 
of K, each time starting with exactly the same initial field 
configuration. It is easy to see that increasing k will lead 
to faster rolling of 9 inside the bubbles. From this one 
may have the impression that larger values of k may lead 
to smaller number of defects as 9 will be close to zero 
in all the bubbles. As mentioned in ref. [^, this indeed 
does happen for values of k which are so large that 9 in 
the bubbles settles down to zero before the oscillations 
can take place. For our choices of parameters in Eqn.(3), 
this happens for values of k larger than 0.03, see ref. [p|. 

However, as long as k is smaller than this value, an 
increase in its value leads to increase in the defect pro- 
duction. It is easy to understand why this happens. In- 
creasing K leads to increase in the potential energy, as 
well as the gradient energy of the field configuration in 
the coalesced region. For example, the region where 9 
takes value tt has much larger energy with larger value of 
K. There is, thus, more energy available for wall to de- 
cay in vortex-antivortex pairs leading to larger number 
of defects. 

The case of twelve bubbles which we described above 
corresponded to the value of k equal to 0.015. We re- 
peated that simulation (with exactly same initial field 
configuration) for two other values of k, one with k equal 
to 0.01 and the other with k equal to 0.02. For the case 
with K equal to 0.01 we got a total of 12 vortices and an- 
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tivortices giving the average number of defects per bub- 
ble to be 1.0. This should be compared with the defect 
density per bubble = 1.5 for k = 0.015 case. From the 
above discussion, this decrease in defect production for 
smaller value of k is expected. For k — 0.02 we found a 
total of 22 vortices and antivortices giving the density of 
defects per bubble to be equal to 1.8, clearly a significant 
increase. In this series of simulations for different val- 
ues of K we found that some of the two bubble collisions 
lead to vortex-antivortex pairs only for k = 0.02. (Note 
that initial bubble configurations have been chosen to be 
identical for all the three cases.) Which implies that (p 
oscillations were large enough to flip $ only for this larger 
value of K. 

There is one important feature we found which needs 
to be emphasized. In the case of k = 0.01 (for twelve 
bubble case), we found that out of the total of 12 vor- 
tices, eight vortices can be attributed to be forming from 
the Kibble Mechanism and four from the pair production 
process. However, when the simulation was repeated for 
K — 0.0 (again starting with exactly same initial field con- 
figuration) then only four vortices formed via the Kibble 
mechanism. The reason that the vortices corresponding 
to the Kibble mechanism are fewer in this case is that 
these corresponded to collisions of more than three bub- 
bles. Altogether these bubbles had zero net winding but 
taken in groups of three bubbles at a time, they had 
winding and anti- winding (using the geodesic rule). In 
the case of k — 0.0 this vortex and antivortex immedi- 
ately annihilate each other (or in some sense their wind- 
ings annihilate each other even before they could form). 
However, for non-zero n, this vortex and antivortex get 
separated in order to minimize the energy as 9 in between 
the vortex and antivortex takes the value zero which is 
the absolute minimum of the effective potential. Thus 
explicit symmetry breaking, in addition to producing de- 
fects via wall decaying into pairs, also can enhance the 
defect production via the Kibble mechanism in certain 
situations. 

To summarize the results of this section, we found that 
increasing the nucleation rate decreases the number den- 
sity of vortices per bubble which is due to decrease in 
average separation between the bubbles. Smaller sepa- 
ration between bubbles leads to a suppression in os- 
cillations and consequently decreases the defect produc- 
tion. Increase in the magnitude of the explicit symmetry 
breaking k, increases the number density of vortices per 
bubble. This is attributed to enhancement in cj) oscil- 
lation coming from the increase in potential energy and 
gradient energy of wall configuration in the coalesced re- 
gion due to a larger k. 

V. EFFECT OF DAMPING ON THE 
NEW MECHANISM 

So far our discussion (as well as the discussion in ref. 
[H) was for the case with no dissipation present. Since 



the creation of vortex-antivortex pairs via this mecha- 
nism crucially depends on oscillations (and on gradi- 
ent which comes from 9 oscillations), it is natural to ex- 
pect that the presence of damping can crucially affect the 
effectiveness of this mechanism. It is thus important to 
understand the precise manner in which damping affects 
this mechanism. After all, damping is naturally present 
in most cases of interest, such as the early Universe, liq- 
uid crystal systems and even quark-gluon plasma at finite 
temperature (for the case of Skyrmions) . 

In this section we will present results of simulations of 
two bubble collisions in presence of dissipation and study 
how the defect production is affected. What we expect to 
find is that the presence of damping should suppress the 
defect production via this mechanism. This is because 
damping will lead to suppression in the amplitudes of 
successive oscillations of (f) in the wall so that it will be- 
come difficult for $ to flip. Consequently, production of 
vortex-antivortex pairs will be suppressed. This is what 
we find in our simulations, which we describe in the fol- 
lowing. 

We have studied the effect of damping on the creation 
of vortices in two bubble collision by introducing a damp- 
ing term rjcj) in the equation of motions for evolving the 
field configuration. Lattice size is taken to be 160 x 128, 
with the same lattice spacing as before. We choose the 
initial field configuration of the two bubbles to be the 
one which leads to formation of three vortex-antivortex 
pairs when damping is absent. Value of k is taken to be 
equal to 0.015. We then repeat the simulation for this 
two bubble collision starting with exactly the same initial 
configuration, but now in the presence of damping. We 
study the effect of magnitude of damping by changing the 
damping coefficient 77. In the presence of damping, bub- 
ble wall velocity is smaller so the bubble collision is less 
energetic to begin with. Further, as the oscillating in 
the coalesced region looses energy also due to damping, 
successive oscillations are strongly suppressed. This leads 
to suppression in the number of pairs getting formed. 

For small values of r/ (in the range 0.0 to 0.25) we still 
get three vortex-antivortex pairs. But for the values of 
r] larger than 0.25, (j) oscillations are strongly suppressed 
and only one pair is formed. For 77 = 1.0, even the one 
pair formed corresponds to a vortex and an antivortex 
which are very close to the wall of the coalesced bubbles 
and are not well formed, cj) oscillations are very strongly 
suppressed in this case and 6 almost settles down to zero 
without any oscillation about 9 = 0.0. Figs. 5a and 5b 
show surface plots of —0 of the coalesced region of the 
two bubbles, at the same stage, t — 84.85, for the two 
cases, rj — 0.1 and 77 — 0.75 respectively. Starting 9 in 
both the bubbles was such that 9 in between the bubbles 
interpolates through 9 ~ tt. As we mentioned above, the 
initial field configurations of the two bubbles were iden- 
tical for both of these cases. We see that there are two, 
well formed, vortex-antivortex pairs present in Fig. 5a, for 
rj = 0.1 case. In contrast, at the same stage, there is only 
one pair present in Fig. 5b for 77 = 0.75 case. For 77 = 0.1 
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case, one more pair forms later, making total number of 
pairs formed equal to 3 for this case. However, by that 
time, the only pair formed for the rj = 0.75 case escapes 
out of the walls of the coalesced bubble (due to 9 being 
zero between the pair, and ^ = tt in directions towards 
the walls), hence we do not show plots at that stage. 




FIG. 5. (a) Surface plot of —0 for rj = 0.1 of the coa- 
lesced region of two bubbles, showing profiles of two vor- 
tex-antivortex pairs. One more pair forms later in this case, 
(b) Surface plot of —0 for rj = 0.75 showing the profile of the 
single pair formed in this case. Later, this vortex and the an- 
tivortex escape out of the coalesced bubble. Both these plots 
are shown at the same time stage, t = 84.85, starting with 
identical initial field configurations of the two bubbles. 

We emphasize here that whatever be the value of the 
damping coefficient, there is always at least one vortex- 
antivortex pair formed with favorable initial configura- 
tion of two bubbles (i.e. 9 distribution) via the new 
mechanism. This happens because the potential energy 
contained in the coalesced region will make the flipping 
of $ (and hence the pair production) favorable compared 
to rolling down of the two ends of the arc in Fig. 1 as that 
leads to constantly increasing gradient energy. We also 
mention that irrespective of damping (or even the or- 
der of phase transition) one pair should always form, for 
appropriate values of 9 in the two bubbles, via the mech- 
anism discussed in ref. [T^ ] as oscillations play no role 
there. However, we have checked that rotation of 9 in the 
walls is so slow (due to the presence of (j)^ in the explicit 
symmetry breaking term) that one will expect the pair 
formation via the mechanism of ref. jTst only after a very 



long time. On the other hand, pair formation due to 
flipping happens in a much shorter time scale, that of a 
single oscillation in the wall. 

Another important point which should be mentioned 
here is that, in the presence of strong damping, no pair 
may get created if k is too small. In the presence of strong 
damping one can neglect the kinetic energy of the bub- 
ble wall. So what remains is the contribution from the 
explicit symmetry breaking term, which must lead to the 
flipping of by forcing it to go through zero. However, 
below a certain critical value of k it may not happen, as 
the 9 gradient developed in the coalesced portion of the 
bubbles may form a stable domain wall. This will be like 
the axionic domain wall ]l^, which can decay only via 
quantum tunneling (say, at zero temperature). 

This is an appropriate point to discuss an important 
issue about the applicability of the new mechanism for 
different types of defects. There is one class of defects 
where the applicability of this mechanism may hold only 
in some special situations, and not generically. These cor- 
respond to vacuum manifolds for which the opposite ori- 
entations of the order parameter field are identified. Ne- 
matic liquid crystals happens to belong to this category, 
with the vacuum manifold being RP^. In such cases, flip- 
ping of the order parameter field does not change its con- 
figuration. Thus the argument given above for the pair 
production can not be directly applied here. However, it 
is possible to argue that in the presence of explicit sym- 
metry breaking this mechanism should still be applica- 
ble, especially when the system is dissipative. Consider, 
for example, an order parameter configuration varying 
around a point P in the physical space smoothly such 
that the value of the order parameter at P is energeti- 
cally most unfavorable (due to explicit symmetry break- 
ing). For the case described by Eqn.(l) it means 6 having 
value TT at point P, and varying from a value less than 
TT to a value larger than tt as we cross P in the physical 
space. 

For nematic liquid crystals, say, in the presence of ex- 
ternal electric field along the x axis, it means that the 
director (order parameter) will lie in the y-z plane at the 
point P, and will deviate from the y-z plane in oppo- 
site directions as one passes through P in the physical 
space. As the order parameter field rolls down to the 
direction of the true vacuum, a region of large gradient 
energy will arise near P. As we explained above, for suf- 
ficiently small explicit symmetry breaking, the gradient 
energy may not become too large, and may result in the 
formation of a domain wall in that region. However, for 
larger values of explicit symmetry breaking (e.g., stronger 
electric field for the liquid crystal case) , the gradient en- 
ergy will keep rising, eventually forcing the field to go 
to zero magnitude in the region near P. The only way 
to decrease the energy of this configuration is by creat- 
ing a defect-antidefect pair. For extremely damped sit- 
uations it may happen gradually as the director reaches 
zero magnitude at P and then slowly grows in magnitude 
but now with an orientation along the x axis, thereby 
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resulting in the creation of vortex-antivortex pair (for 2- 
dimensions, and a loop of string for 3-dimensions). For 
smaller damping, or for that matter, in the absence of 
damping, the order parameter field may keep oscillating 
at P through zero, without changing its configuration in 
passing through zero. During the passage of the field 
through zero, any small fluctuation can re-orient the di- 
rector along the x axis, resulting in the pair production. 
Or, the field may eventually settle at zero (either due to 
the damping term, or just due to loss of energy in succes- 
sive oscillations in other excitations of the system), and 
then roll down to a direction along the x axis. This will 
then again result in the pair production via this mecha- 
nism. In the absence of any explicit symmetry breaking, 
this mechanism does not seem to be applicable for these 
types of vacuum manifolds. 

VI. CONCLUSIONS 



not even a single defect is produced via the Kibble mech- 
anism, while 104 vortices and antivortices are produced 
via the new mechanism. As the distribution of these 
defects is of very different nature than the one expected 
from the Kibble mechanism, one may expect qualitatively 
new features for systems with explicit symmetry break- 
ing. In our 2+1 dimensional study we find a sequence 
of vortex-antivortex pairs being produced. Simple argu- 
ments show that for 3-1-1 dimensions one will get a system 
of concentric string loops, (see ref. [|j for example). This 
may then suppress formation of large strings compared 
to small loops. It will be of great interest if such a pre- 
diction can be experimentally verified in some physical 
system. As mentioned in ref. formation of strings in 
nematic liquid crystals in the presence of external electric 
or magnetic fields may be ideal for checking this mech- 
anism. (Though in that case, presence of damping may 
suppress string formation via this new mechanism.) 



We conclude by emphasizing the essential aspects of 
our results. We have carried out numerical simulations 
of a first order phase transition by random nucleation 
of bubbles, in the presence of small explicit symmetry 
breaking, and have studied the production of vortices and 
antivortices. We estimate the net number density (num- 
ber of defects per bubble) of vortices produced, which 
includes vortices formed due to the Kibble mechanism as 
well as those produced via the pair production mecha- 
nism. We also study the dependence of this defect num- 
ber density on parameters such as the magnitude of the 
explicit symmetry breaking term as well as on the nu- 
cleation rate. Defect production increases with larger 
magnitude of explicit symmetry breaking due to larger 
potential energy in the coalesced region. Nucleation rate 
affects defect density due to the fact that a larger nu- 
cleation rate implies smaller average bubble separation, 
which in turn leads to less kinetic energy for the bub- 
ble walls before bubbles collide. Oscillations of are less 
prominent for less energetic walls leading to smaller num- 
ber of defects for larger nucleation rates. We have also 
studied the effects of presence of damping on this pair 
production mechanism by studying pair production in 
two bubble collisions. As expected, we find the damping 
suppresses oscillations of and hence leads to smaller 
number of defects, though one pair is always produced 
(for suitable values of 9 distribution in the bubbles, and 
for K not too small). 

In all these cases we find that the number of defects 
produced due to effects of explicit symmetry breaking 
term (either via direct pair production mechanism, or by 
separating vortex-antivortex pairs produced in collisions 
of more than three bubbles via the Kibble mechanism) is 
much larger than what one would expect from the Kibble 
mechanism. This relative enhancement is much larger 
when we include the effect of bias in distribution in 
bubbles, resulting from explicit symmetry breaking. In 
this case we find, in a simulation with 63 bubbles, that 
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